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The thermal evolution of the cosmic gas decoupled from that of the cosmic microwave background 
(CMB) at a redshift z ~ 200. Afterwards and before the first stars had formed, the cosmic neutral 
hydrogen absorbed the CMB flux at its resonant 21cm spin-flip transition. We calculate the evolu- 
tion of the spin temperature for this transition and the resulting anisotropics that are imprinted on 
the CMB sky due to linear density fluctuations during this epoch. These anisotropics at an observed 
wavelength of 10.56 {(1 + z)/50} meters, contain an amount of information that is orders of magni- 
tude larger than any other cosmological probe. Their detection, although challenging, could tightly 
constrain any possible running of the spectral index from inflation (as suggested by WMAP), small 
deviations from Gaussianity, or any significant contribution from neutrinos or warm dark matter to 
the cosmic mass budget. 

PACS numbers: 98.80.-k, 98.65.-r, 98.70.Vc, 95.30.Jx 

Introduction. The small residual fraction of free elec- 
trons after cosmological recombination coupled the tem- 
perature of the cosmic gas to that of the cosmic mi- 
crowave background (CMB) down to a redshift, z ~ 200 
[1] . Subsequently, the gas temperature dropped adiabat- 
ically as Tgas oc (1 -I- z)^ below the CMB temperature 
T-y (X (1 -I- z). The gas heated up again after being ex- 
posed to the photo-ionizing ultraviolet light emitted by 
the first stars during the reionization epoch at z ^ 20 (see 
review in [2]). Prior to the formation of the first stars, the 
cosmic neutral hydrogen must have resonantly absorbed 
the CMB flux through its spin-flip 21cm transition [3-6]. 
The linear density fluctuations at that time should have 
imprinted anisotropics on the CMB sky at an observed 
wavelength of 21.12[(1 -I- 2:)/100] meters. In this Letter, 
we calculate the power spectrum of these anisotropics and 
assess the significance of their potential detection. A di- 
rect measurement of the amplitude of inhomogeneities on 
small spatial scales would constrain any possible tilt or 
running of the spectral index of the power spectrum as 
recently suggested by Wilkinson Microwave Anisotropy 
Probe ( WMAP) [7] , or any suppression of power on small 
scales due to a warm dark matter component in the cos- 
mic mass budget [8]. 

Spin temperature history. We start by calculating the 
history of the spin temperature, Tg, defined through 
the ratio between the number densities of hydrogen 
atoms in the excited and ground state levels, ni/no = 
(51/50) exp {—T^/Ts} , where subscripts 1 and corre- 
spond to the excited and ground state levels of the 21cm 
transition, (gi/go) = 3 is the ratio of the spin degener- 
acy factors of the levels, nn = (?^o + '^1) c>c (1 -|- z)^ is the 
total hydrogen density, and T^, = 0.068K is the tempera- 
ture corresponding to the energy difference between the 
levels. The time evolution of the density of atoms in the 
ground state is given by. 
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where a{t) = (1 + z)^^ is the cosmic scale factor, A's 
and S's are the Einstein rate coefficients, C's are the 
coUisional rate coefficients, and is the blackbody in- 
tensity in the Rayleigh- Jeans tail of the CMB, namely 
= 2kBT^/\^ with A = 21 cm [9]. The 0^1 transi- 
tion rates can be related to the 1 — + transition rates 
by the requirement that in thermal equilibrium with 
Ts — T-y = Tgas, the right-hand-side of Eq. (1) should 
vanish with the coUisional terms balancing each other 
separately from the radiative terms. The Einstein coef- 
ficients are ^10 = 2.85 x IQ-^^ s"\ ^ {X^/2hc)Aio 
and i3oi = (51/50)^10 [3,9]. The coUisional de-excitation 
rates can be written as Cio = |k(1 — 0)nH, where k(1 — 0) 
is tabulated as a function of Tgas [10]. 

Equation (1) can be simplified to the form. 
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where T = no/nn, H w i/o\/^m(l + z)^/^ is the Hubble 
parameter at high redshifts (with a present-day value of 
Ho), and flm is the density parameter of matter. The 
upper panel of Fig. 1 shows the results of integrating 
Eq. (2). Both the spin temperature and the kinetic tem- 
perature of the gas track the CMB temperature down to 
z ~ 200. Collisions are efficient at coupling T, and Tgas 
down to z 70 and so the spin temperature follows the 
kinetic temperature around that redshift. At much lower 
redshifts, the Hubble expansion makes the collision rate 
subdominant relative the radiative coupling rate to the 
CMB, and so Tg tracks Ty again. Consequently, there is 
a redshift window between 30 ^ z < 200, during which 
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the cosmic hydrogen absorbs the CMB flux at its reso- 
nant 21cm transition. Coincidentally, this redshift inter- 
val precedes the appearance of collapsed objects [2] and 
so its signatures are not contaminated by nonlinear den- 
sity structures or by radiative or hydrodynamic feedback 
effects from stars and quasars. 
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FIG. 1. Upper panel: Evolution of the gas, CMB and spin 
temperatures with redshift [4]. Lower panel: dTb/dSn as func- 
tion of redshift. The separate contributions from fluctuations 
in the density and the spin temperature are depicted. We also 
show dTt/dSffCi oc dTi/dSn x Sh, with an arbitrary normal- 
ization. Throughout this Letter, we assume the standard set 
of cosmological parameters for a universe dominated by cold 
dark matter and a cosmological constant (ACDM) [6]. 

During the period when the spin temperature is 
smaller than the CMB temperature, neutral hydrogen 
atoms absorb CMB photons. The resonant 21cm ab- 
sorption reduces the brightness temperature of the CMB 
by, 

Tb = T(r,-T^)/(l + z), (3) 
where the optical depth for resonant 21cm absorption is, 



32nkBTsH{z) 
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Small inhomogeneities in the hydrogen density Sh = 
(rifj — n}i)/nii result in fluctuations of the 21cm absorp- 
tion through two separate effects. An excess of neutral 
hydrogen directly increases the optical depth and also al- 
ters the evolution of the spin temperature. We can write 
an equation for the resulting evolution of T fluctuations. 



dST 

dz 



= [Hil + z)]-^ {[Cio + C'oi + (Boi + Bio)I,]6T 



+ [CoiT-Cw{1-T)]Sh}, 
leading to spin temperature fluctuations, 
6Ts 1 ST 



ln[3T/(l-T)]T(l-T)' 
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The resulting brightness temperature fluctuations can be 
related to the derivative, 



dSn 
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through 6Tb = {dTb/d6H)SH- We include all fluctua- 
tions caused by Sh except for the variation in due 
to fluctuations in Tgas which is very small [10]. Figure 
1 shows dTb/dSu as a function of redshift, including the 
two contributions to dTb/dSn, one originating directly 
from density fluctuations and the second from the associ- 
ated changes in the spin temperature [4] . Both contribu- 
tions have the same sign, because an increase in density 
raises the collision rate and lowers the spin temperature 
and so it allows Ts to better track Tgas- Since Sh grows 
with time as Sh oc a, the signal peaks at z 50, a slightly 
lower redshift than the peak of dTb/dSn- 

Next we calculate the angular power spectrum of the 
brightness temperature on the sky, resulting from density 
perturbations with a power spectrum Ps{k), 



(<5H(ki)<5H(k2)) = (27r)^^^(ki +k2)P5(fc 
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where 5//(k) is the Fourier tansform of the hydrogen den- 
sity field, k is the comoving wavevector, and (• • •) denotes 
an ensemble average (following the formalism described 
in [6]). The 21cm brightness temperature observed at a 
frequency v corresponding to a distance r along the line 
of sight, is given by 

STb{n,i.) = J drWAr) ^ShM, (9) 

where n denotes the direction of observation, W^{r) is 
a narrow function of r that peaks at the distance cor- 
responding to V. The details of this function depend 
on the characteristics of the experiment. The brightness 
fluctuations in equation (9) can be expanded in spherical 
harmonics with expansion coefficients a/„i(i^). The angu- 
lar power spectrum of map Ci{v) — (|a;m(t^)p) can be 
expressed in terms of the 3D power spectrum of fluctua- 
tions in the density Ps{k), 
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Our calculation ignores inhomogeneities in the hydrogen 
ionization fraction, since they freeze at the earlier recom- 
bination epoch (z ~ 10"^) and so their amplitude is more 
than an order of magnitude smaller than at z < 100. 
The peculiar velocity and gravitational potential pertur- 
bations induce redshift distortion effects that are of order 
~ (H/ck) and ^ (H/ck)^ smaller than Sh for the high-? 
modes of interest here. These effects are expected to be 
washed-out within realistically broad band filters. 
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FIG. 2. Angular power spectrum of 21cm anisotropies 
on the sky at various redshifts. From top to bottom, 
z = 55,40,80,30, 120,25,170. 



Figure 2 shows the angular power spectrum at var- 
ious redshifts. The ability to probe the small scale 
power of density fluctuations is only limited by the Jeans 
scale, below which the dark matter inhomogeneities are 
washed out by the finite pressure of the gas. Interestingly, 
the cosmological Jeans mass reaches its minimum value, 
^ 3 X 10* Mq, within the redshift interval of interest here 
[2]. During the epoch of rcionization, photoionization 
heating raises the Jeans mass by several orders of mag- 
nitude and and broadens spectral features, thus limiting 
the ability of other probes of the intergalactic medium, 
such as the Lya forest, from accessing the same very low 
mass scales. The 21cm tomography has the additional 
advantage of probing the majority of the cosmic gas, in- 
stead of the trace amount 10^^) of neutral hydrogen 
probed by the Lya forest after rcionization. Similarly to 
the primary CMB anisotropies, the 21cm signal is simply 
shaped by gravity, adiabatic cosmic expansion, and well- 
known atomic physics, and is not contaminated by com- 
plex astrophysical processes that affect the intergalactic 
medium at z < 30. 

The small scale power spectrum. In most models of 
inflation, the evolution of the Hubble parameter during 
inflation leads to departures from a scale-invariant spec- 
trum that are of order l/A'cfoid with iVofoid ~ 60 being 
the number of e-folds between the time when the scale of 
our horizon was of order the horizon during inflation and 
the end of inflation [14]. Recent WMAP data combined 
with other measures of the power on smaller scales, sug- 
gests that the power spectrum changes with scale much 
faster than inflation would have predicted [7], although 
this result is still somewhat controversial. Independent 
hints that the standard ACDM model may have too much 
power on galactic scales have inspired several proposals 
for suppressing the power on small scales. Examples in- 
clude the possibility that the dark matter is warm and it 



decoupled while being relativistic so that its free stream- 
ing erased small-scale power [8], or direct modiflcations 
of inflation that produce a cut-off in the power on small 
scales [15]. An unavoidable coUisionless component of the 
cosmic mass budget beyond CDM, is provided by mas- 
sive neutrinos (see [16] for a review). Particle physics 
experiments established the mass splittings among dif- 
ferent species which translate into a lower limit on the 
fraction of the dark matter accounted for by neutrinos of 
fi, > 0.3%, while current constraints based on galaxies 
as tracers of the small scale power imply < 12% [17]. 

In Fig. 3 we show the 21cm power spectrum for various 
models that differ in their level of small scale power. It 
is clear that a precise measurement of the 21cm power 
spectrum will dramatically improve current constraints 
on alternatives to the standard ACDM spectrum. 




FIG. 3. Upper panel: Power spectrum of 21cm anisotropies 
at z = 55 for a ACDM scale-invariant power spectrum, 
a model with n = 0.98, a model with n = 0.98 and 
Or = In P/dln fc^) = —0.07, a model of warm dark mat- 
ter particles with a mass of 1 keV, and a model in which 
fv ~ 10% of the matter density is in three species of massive 
neutrinos with a mass of 0.4 eV each. Lower panel: Ratios 
between the different power spectra and the scale-invariant 
spectrum. 

Unprecedented information. The 21cm signal con- 
tains a wealth of information about the initial fluctua- 
tions. A full sky map at a single photon frequency mea- 
sured up to Zmax, can probe the power spectrum up to 
^max ^ (^max/10^)Mpc~^. Such a map contains l^^^y, in- 
dependent samples. By shifting the photon frequency, 
one may obtain many independent measurements of the 
power. When measuring a mode I, which corresponds 
to a wavenumber k ~ l/r, two maps at different photon 
frequencies will be independent if they are separated in 
radial distance by 1/k. Thus, an experiment that covers 
a spatial range Ar can probe a total of kAr ~ lAr/r in- 
dependent maps. An experiment that detects the 21cm 
signal over a range Ai/ centered on a frequency v, is sensi- 
tive to Ar/r ~ 0.5(Aj//j/)(1 -|- and so it measures 
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atotalof TVaicm - 3xlO^^{U^^/10^f{Au/u){z/100)-'^/'^ 
independent samples. 

This detection capability cannot be reproduced even 
remotely by other techniques. For example, the pri- 
mary CMB anisotropics are damped on small scales 
(through the so-called Silk damping), and probe only 
modes with I < 3000 {k < 0.2 Mpc"^). The total num- 
ber of modes available in the full sky is iVcmb = 2^^ax ~ 
2 X 10^(Zmax/3000)^, including both temperature and po- 
larization information. 

Detectability of signal. The sensitivity of an experi- 
ment depends strongly on its particular design, involving 
the number and distribution of the antennae for an in- 
terferometer. Crudely speaking, the uncertainty in the 
measurement of [l{l + l)C;/27r]^/^ is dominated by noise, 
N,y, which is controlled by the sky brightness at the 
observed frequency u [6], 



7V^~ 0.4mK 



"35" 



5000 



0.016 



fct 



5 X 10^ Jy Sr ^/ \ 35 J \ /max / \^cover 



1 \ 1/2 / A 

1 year \ / 
to 



(11) 



where Zmin is the minimum observable I as determined 

by the field of view of the instruments, ^niax is the max- 
imum observable I as determined by the maximum sep- 
aration of the antennae, /cover is the fraction of the ar- 
ray area thats is covered by telescopes, to is the obser- 
vation time and Aiy is the frequency range over which 
the signal can be detected. The numbers adopted above 
are appropriate for the inner core of the LOFAR ar- 
ray {http://www.lofar.org), planned for initial opera- 
tion in 2006. The predicted signal is ~ ImK, and 
so a year of integration or an increase in the cover- 
ing fraction are required to observe it with LOFAR. 
Other experiments whose goal is to detect 21cm fluc- 
tuations include SKA {http://www.skatelescope.org) and 
PAST {http://astrophysics.phys. emu. edu/^jbp). The 
main challenge in detecting the predicted signal involves 
its appearance at low frequencies where the sky noise 
is high. Proposed space-based instruments [18] avoid 
the terrestrial radio noise and the increasing atmospheric 
opacity at ^ 20 MHz (corresponding to z > 70). 

Final comments. The 21cm absorption is replaced by 
21cm emission from neutral hydrogen as soon as the in- 
tergalactic medium is heated above the CMB tempera- 
ture during the epoch of reionization [13]. Once most of 
the cosmic hydrogen is reionized at Zreiom the 21cm signal 
is diminished. The optical depth for free-free absorption 
after reionization, ~ 0.1[(l-|-2;reion)/20]^/^, modifies only 
slightly the expected 21cm anisotropics. Gravitational 
lensing should modify the power spectrum [12] at high I, 
but can be separated as in standard CMB studies (see [19] 
and references therein). The 21cm signal should be sim- 
pler to clean as it includes the same lensing foreground 
in independent maps obtained at different frequencies. 



The large number of independent modes probed by the 
21cm signal would provide a measure of non-Gaussian 

— 1/2 

deviations to a level of ~ iV2icm > constituting a test of 
the inflationary origin of the primordial inhomogencities 
which are expected to possess deviations > 10~^ [20]. 
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